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We explore quantum search from the geometric viewpoint of a complex projective space CP, a 
space of rays. First, we show that the optimal quantum search can be geometrically identified with 
the shortest path along the geodesic joining a target state, an element of the computational basis, 
and such an initial state as overlaps equally, up to phases, with all the elements of the computational 
basis. Second, we calculate the entanglement through the algorithm for any number of qubits n 
as the minimum Fubini-Study distance to the submanifold formed by separable states in Segre 
embedding, and find that entanglement is used almost maximally for large n. The computational 
time seems to be optimized by the dynamics as the geodesic, running across entangled states away 
from the submanifold of separable states, rather than the amount of entanglement itself. 

PACS numbers: 03.67.Lx, 03.65.-w, 89.70.+C 



I. INTRODUCTION 

Quantum computers would be more powerful than 
their classical counterparts [Q,^] . Suppose an oracle func- 
tion f(x) with x € {0, 1}" is given such that f(w) = 1 for 
an unknown single item w out of N(:= 2") and f(x) = 
for x ^ w. Our purpose is to find the "target" w with the 
smallest possible number of the oracle evaluations, called 
the query-complexity. As is often the case with computer 
science, the worst case of query-complexity is concerned 
here. If we try with a classical computer, it is readily 
found that we need TV queries in the worst case. On the 
other hand, we can obtain w with the success probability 
almost 1 in only 0(\HV) queries, regardless of w (i.e. for 
the evaluation in not only the worst case but also the av- 
erage case), by Grover's quantum search algorithm 
Furthermore, Zalka [^) proved that Grover's algorithm is 
exactly, and not only asymptotically, optimal for query- 
complexity if quantum computation consists only of uni- 
tary transformations and the final measurement. 

Grover's algorithm in rc-qubits (2" = N states) case 
is constructed as follows. We first introduce an initial 
"average" state \a) := Y^,x=o \ x ) wnere \ x ) ( x = 
0, . . . , N— 1) forms the orthonormal computational ba- 
sis. Writing the overlap between the average \a) and the 
target \w) by 6 as 
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we have Grover's algorithm: 
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in the orthonormal basis of \r)(:= -^=^2 x _ iw \x)) and 

\w), where 1 denotes the 2nx2n identity matrix. Note 
that by constructed from alternate inversions for the av- 
erage la) and the target \w) (i.e. I a and I w ), the kernel 
G in (ph becomes a real 2-dimensional rotation. We find 
in (|H) that the target \w) is obtained with the success 
probability 1 when {k+\)9 = tt/2, i.e. k ~ f \/F by (Q) 
in case oi N ^> 1. Because one query is used for every 
I w (i.e., I w \x) = (-l) f W\x)) of G in (|), we can thus 
identify the query-complexity with k ~ 0(y/~N). 

Now our motivation is based on the two points below: 
first, while it is quite straightforward to verify Grover's 
algorithm 11,0] and Zalka's algebraic proof of its opti- 
mality 0], it has yet to be understood from the geo- 
metric aspects why Grover's algorithm works efficiently. 
Second, though it is often said that entanglement is use- 
ful to enhance the quantum information processing, this 
remains obscure in theory as well as in the Nuclear 
Magnetic Resonance (NMR) experiments (9| after, in par- 
ticular, Lloyd's proposition of "quantum search without 
entanglement" [10,0. Thus, in this paper, we char- 
acterize quantum search from the geometric viewpoint, 
which might shed light on the general strategy for con- 
structing efficient quantum algorithms, and discuss how 
entanglement gives quantum computation its power. 

The rest of the paper is organized as follows. In Sec. |J, 
after we briefly review geometric aspects of quantum me- 
chanics such as the complex projective space CP and 
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Fubini- Study metric on it, we show that Grovcr's algo- 
rithm corresponds to a geodesic of CP. In Sec. [II wc dis- 
cuss entanglement, which can be considered as the min- 
imum Fubini-Study distance to the submanifold formed 
by separable states in CP. Entanglement in Grover's al- 
gorithm is calculated for the general n-qubits case, and 
is found to be used almost maximally when n is large. In 
Sec. IV, we construct optimal quantum searches, includ- 
ing Grover's algorithm, by means of geodesies, and derive 
the geometric necessary and sufficient condition for the 
optimal quantum search. Finally Sec. ^ is devoted to 
conclusions. 



called the inhomogeneous coordinates. 



B. Fubini-Study metric and geodesies in CP 

Now wc introduce Fubini-Study metric, a natural Rie- 
mannian metric in CP _1 . Let \tp(s)) be a normal- 
ized vector drawing a curve C in Ti and \dip(s)) the 
tangent vector along C. Note that the normalization: 
(ip(s)\ip(s)) = 1 implies HB(i/j(s)\dijj(s)) = 0. Under a 
global gauge transformation: \tp) i— > e n \ip) with 7 G M, 
the projection, orthogonal to the Hopf fibers, of \dijj(s)): 



II. GEOMETRIC ASPECTS OF QUANTUM 
MECHANICS 

In this section, we first consider, for preliminaries, the 
pure state space of a quantum mechanical system as a 
complex projective space CP, a space of rays in the as- 
sociated Hilbert space Ji ]l2|-|l8fl . Because we discuss the 
geometric characters of the efficient quantum algorithm 
itself, we can safely restrict to our attentions to the pure 
states. It implies that we never treat general mixed states 
(whole states given by the density matrix) which appear 
in some realistic situations. After that, we show that 
Grover's algorithm is the horizontal lift of a geodesic in 
CP. 



A. Ray and complex projective Hilbert space CP 

Let be a (not necessarily normalized) vector in a 
complex A^-dimensional Hilbert space H(C N ). The phys- 
ical state of the quantum system in H(C ) is given by a 
ray, an equivalence class of vectors up to the overall nor- 
malization and phase. So the ray can be interpreted as 
a line in C N passing through the origin. Note that uni- 
versal quantum computation [?| is defined over rays. 
A set of rays forms the complex projective Hilbert space 
qypN-i w jth the associated projection map n, 
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{W) s.t. |V) =<#>, ceC-{0}}. 
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Suppose \ip) is given by A^-tuplcs of complex amplitudes 
zj (J = 0, ... ,N—1) G C^ — {0} by choosing a basis in 
Ji. According to (|^), the ray n(j?/>)) is represented as 

n(|V}) = ••• >4r-i). ( 6 ) 



such that z'a 



cz 



for all j with c G C-{0}. We find, 



however, this representation (g), called the homogeneous 
coordinate representation in algebraic geometry, is not 
unique. To obtain a unique one, we also utilize, for any 
nonzero Zj (say zq), 







-H- (1 = 1,. 



--0 



jv-i), 



(7) 



|#(sU) := |#(s)) - (^(s)|#( S ))|V(s)) (8) 



is gauge covariant (i.e. \dip±) 



e^\dip ± )). Since 



(dip(s)j_\dip(s)±) is gauge invariant, it can be used to 
define the metric in CP^ -1 , called Fubini-Study metric, 
between two nearby rays TL(\ip(s)}) and H(\ip(s+ds))) as 



ids 2 := (<ty(*)x|#(*)x> 

= (d^( s )\d^(s)) - (Im(lK*)|#(*)» S 



(9) 



By variation of the action J* 2 ds of the line element 
in (||), each extremal gives a geodesic C , which is found 
to be an arc of the great circle lying on some submani- 
fold CP 1 in <CV N - 1 Any lift of the geodesic C 
becomes, by definition, a geodesic in Ji. In particular, a 
horizontal lift of C' , which implies the parallel transport: 
lm( , ip(s)\d'ip(s)} = 0, can be described simply as 



|V>(s)) = cos- l^i) 



sin -\ijj 2 ) 



(10) 



in terms of some orthonormal basis \ipi), i n Thus 
the horizontal geodesic (jlTJ) is just a real 2-dimensional 
rotation on the plane spanned by and \ip2) in H. Fur- 
thermore, according to (pTj|), we can interpret the transi- 
tion probability P as the distance s (G [0,7r]) along the 
geodesic joining \i/)x) and \tp{s)) fl^Hlql; i.e. 



pm S )),\ipi))-.= \msMi 



(ii) 



We also find the geodesic represents possible superposi- 
tions between \ipi) and \ip2)- 



C. Grover's algorithm as a geodesic 



If we take \tpx) = \r), \ip 2 ) = \w) and s = 2(k+\)6 in 
(|l0|), wc readily find Grover's dynamics (||) satisfies the 
equation of a geodesic in (|l0|), in addition evolves surely 
along the shorter arc of the geodesic. This suggests that 
Grover's dynamics corresponds to the shortest path from 
the geometric viewpoint. 
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It is significant to note that the original Grover's al- 
gorithm evolves with discrete k, in other words, it skips 
along the geodesic. The interval of skip becomes shorter 
as N becomes larger, and almost continuous when N is 
sufficiently large. Here we can regard G in (|J) as a one- 
step time evolution, because we are concerned with the 
computational complexity only in terms of the number 
of queries called. This might be called "coarse-graining," 
where the dynamics driven by the detailed physical op- 
erations is reduced to the effective dynamics (i.e. algo- 
rithm) of the query-complexity. 

III. ENTANGLEMENT IN GROVER'S 
ALGORITHM 

In this section, we explore the geometry of CP in more 
details to consider relationships between Grover's algo- 
rithm and entanglement. As mentioned in Sec. [j], it is 
very interesting whether Grover's algorithm (and general 
quantum algorithms) takes advantage of entanglement to 
compute faster. It is shown in || that Grover's algorithm 
both in the ideal pure state case and in the pseudo-pure 
state in NMR case does generate entanglement during the 
computation, by formally tracing out all but one qubit. 
Here we show from the geometric viewpoint that entan- 
glement is used and calculate it explicitly. 

A. Segre embedding and quadric of separable states 

Some of the mysterious features of quantum mechanics, 
e.g. entanglement and so on, appear when we consider 
a composite system. In the bipartite case, by combin- 
ing two systems with Hilbcrt space TC(C m ) and Tt(C m ), 
the combined Hilbert space is taken as the tensor prod- 
uct H(C m ) <8> H(C m ) and the associated space of states 
is CP mm_1 , which has a much larger dimension than 
that of the mere Cartesian product CP" 1 ' 1 x CP™'" 1 
(its dimension is only m+m' — 2) of the two individual 
spaces of states. Thus the mysteries seem to lie in the 
(m — \){m! — 1) relative phases. Here wc consider Scgrc 
embedding p^-|T^] in algebraic geometry, which enables 
products of projective spaces to be embedded into a pro- 
jective space again. Then using the Segre embedding, we 
may characterize entanglement geometrically. 

We first illustrate the idea in 2-qubits case. (Segre em- 
bedding in the general case is given in Appendix [A|) A 
state of a qubit is represented by the homogeneous co- 
ordinates (zq,Z\) €E CP 1 . In particular, the spin- up and 
spin-down basis states |0) and |1) correspond to 

|0)^(1,0), |1)<->(M), (12) 

respectively (precisely speaking, n(|0)) = (1,0) and 
n(|l)) = (0,1)). Then an arbitrary state (zo,Z\) is a 
point on the complex projective line joining IT( 1 0) ) and 
II(|1)), interpreted as a superposition of |0) and |1) with 



the am pli tudes proportional to zq, z\ respectively as seen 
in Scc.fjJ 

We consider a mapping / (Segre embedding); 

/ : CP 1 x CP 1 -> CP 3 (lS) 
((a ,ai), (&oA)) ^ (oo6o)Oo&i,ai&o,ai6i). 

Note that although (ao,ai) = (aao, crai), (&o,&i) = 
(fibo^h) with a,0 G C-{0}, the above / in @ 
map them to the identical point in CP 3 , regardless of 
a, (3. Now we discuss the condition for the image of / : 
/(CP 1 x CP 1 ) to satisfy in CP 3 . By writing down the 
homogeneous coordinates in CP 3 as (zq, Z\, z%, Z3), we 
define a polynomial of degree 2; 

Q := z z 3 - zjz 2 , (14) 

which satisfies Q(ao&Oi ^o^i, ai&o? a-ibi) = 0. On the 
other hand, it is readily checked that arbitrary points 
on Q = are included in /(CP 1 x CP 1 ). Thus we find 

f(CP 1 xCP 1 )={(z Q ,z 1 ,z 2 ,z 3 ) | Q = 0}. (15) 

Because we can transform any non-singular quadric into 
the "normal" quadric form Q = zqz 3 — ziz 2 = by a 
projective transformation g: 

N-l 

Zj&J2 A V Zl U = 0,...,N-1), (16) 

1=0 

with an A x A matrix A := {A 3 i G C} s.t. dctA ^ 
0, we can also identify the non-singular quadric with 
CP^CP 1 . That is the reason why the algebraic subman- 
ifold of separable, or no-entangled, states (CP 1 x CP 1 ) 
forms the quadric Q = in general state space for the 
2-qubits system (CP 3 ), and the states in CP 3 off the 
quadric Q — are entangled states. 




FIG. 1. The geometry of CP 3 around |11), the assumed 
target \w). Note that because three axes of Ci, C2, and (3 rep- 
resent complex numbers, this figure is written in the complex 
dimension 3 (real dimension 6). By extracting the real axes 
of Cii C2, and £3, 1 11} is found to lie on a saddle point of the 
quadric. So is each of the other states of the computational 
basis. 
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B. Usage of entanglement 

Let us examine the geometry of CP 3 in more detail. 
We suppose the case where the target state \w) is |11) 
without loss of generality. To draw the behavior around 
(0,0,0,1)) as in Fig. [I], we introduce the inho- 
mogeneous coordinates: Co = 20/23, Ci = 24/2:3, and 
C2 = 22/23 because of 23 7^ 0. Then the quadric Q = 



in (14) is written as 



Co = C1C2. 



(17) 



In particular, all the states orthogonal to |11), including 
1 00), |01), 1 10), (|01) + |10))/V2, etc., are located in points 
at infinity in Fig. |l|. The evolution of Grover's algorithm 
in (El) is given by 



cos(A;+^)0 cos(fc+i)# cos(fc+| 
(u,u,u, 1), 



1 



- ? =^, sin(fc+|)# 
(18) 



or in terms of inhomogeneous coordinates 



Co = Ci = C2 



(19) 



where u := cot(fc + i)(9/\/3 changes from 1 to 0. We 
find that Grover's algorithm starts from the average state 
|a) (Co = Ci = C2 = 1) on the quadric, evolves away from 
the quadric along a (geodesic) line < Co = Ci = C2 < 1, 
and finally reaches the target \w) at the origin on it. 
Hence Grover's algorithm uses the entanglement in 2- 
qubits case. 

Now we treat the general n-qubits case. Recalling (Q), 
we now represent \ip(k)) by the homogeneous coordinates 
in CV N - X (N = 2 n ); 



(k) 



cos(fc+i)6» 



z w (k)=sm{k+l)6. (20) 



We discuss whether the states of Grover's evolution in 
CP 2 ~ 1 are included in the algebraic submanifold of 
the completely separable states of CP 1 x ••• x CP 1 (=: 
(C'P 1 ) xn ). As a first step, we consider the condition that 
© arc included in CP 2 " 1-1 x CP 1 . This is just the nec- 
essary condition for the reduction CP 2 _1 — > (CP 1 ) XTt 
and, according to (A3) in Appendix |A| (m = 2"~ 1 — 1 and 
ml = 1), is given by 



cos(k+\)6sm(k+\)6 cos 2 (fc+±)6> 



N 



(21) 



From (pl|), we have two cases: 

(i) If cos(fc + 7^ 0, the condition ( pl| ) becomes 
tan(fc + i)0 = 1/y/N—l. The solutions are given 
as (fc+|)0 = |, § + 7T (mod 2tt) by use of (@). 



(ii) If cos(fc+±)6» = 0, it means {k + \)6 = 2 
(mod 2tt). These are also the solutions of (pi|). 



3tt 
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The solutions (i) and (ii) are also sufficient, i.e. com- 
pletely separable to (CP 1 ) xn , and indeed correspond 
to the average and target state respectively. For the 
states in ( |20| ) with other k, we cannot reduce them into 
(CP 1 ) xn and thus they are entangled states. In brief, 
though the initial (average) state and the target state 
are separable, the intermediate states through which the 
system evolves arc entangled. 



C. Calculation of entanglement 

For the Grover's evolution \ip), let us calculate the 
amount of entanglement E. Entanglement E in our 
pure state space is naturally considered the minimum 
Fubini-Study distance s to the submanifold formed by 
completely separable states (CP 1 ) xn in CP^ -1 , i.e. 



W»:= min s(|V>>, |<£» 

{|0)|(CT>1)><»} 







2arccos / max P(\ip),\6)). 
{|*>|(Cpi)*»} 



(22) 



Because Fubini-Study distance in (|22j) can be reduced 
from Burcs distance with the parallel transport connec- 
tion in the case of pure states |l{| , it satisfies the require- 
ments for a good measure of entanglement: i.e. (i) zero 
for any separable state; (ii) invariant under local unitary 
transformations; (iii) non increasing expectation value 
under local operations, such as classical communication 
and subselection, as Vedral et al. suggested in |2C|| . 

It should be remarked that in the case of a bipartite 
(2-qubits) pure state system, the partial entropy (von 
Neumann entropy of the reduced density matrix associ- 
ated with one of the parties) is widely supposed to be a 
good measure of entanglement pl| . However, we apply 
( p2| ) as the geometric entanglement measure, because (i) 
the partial entropy has no apparent geometric meaning in 
CP; and (ii) an extension to the multipartite (n-qubits) 
case is non trivial p2| . As a comparison, we calculate, 
in Appendix the entanglement by the partial entropy 
in 2-qubits case. We find our measure of entanglement 
( [22| ) almost corresponds to "concurrence" [|3| so as to be 
consistent with the calculation using the partial entropy. 

Let us first discuss the 2-qubits case, and then proceed 
to the general n-qubits case. To calculate the entangle- 
ment of the Grover's state \tp) (in (|l8|)) for 2-qubits case, 
we have to look for the point that gives the minimum 
of s (or maximum of P) in (|2^) on the submanifold of 
the quadric Q = 0. Because this point must lie on the 
plane Ci = C2 (i-e. z\ = 22) as seen in Fig. ^, we can 
parametrize its candidates \<j)) as 



(v 2 ,v,v, 1), 



(23) 
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with v € C such that < \v\ < 1. Thus we consider 
max?(u, v) 

V 

[u(v + l) 2 - u + l)[u(v* + l) 2 - u + 1] 

— IHcLX ^ 

(3u 2 + l)(M 2 + l) 2 
(ur 2 + 2ur cos x + l) 2 + 4u(u — l)r 2 sin 2 % 



max ■ 

(r,x) 



(3u 2 + l)(r 2 + l) 2 



(24) 



where v* denotes the complex conjugate of v and we 
use v = re zx with < r < 1. For fixed r, P in 
is the largest for the phase x — (mod 2ir). This 
is solely because u in ( O ) is a real number. Since 
J^P(u, r, x = 0) = 0, we have, according to 



2 arccos - 



3u 2 



1 



1 



vm 



(25) 



where ijjvj := [u — 1 + y/(u— l) 2 + 4w 2 ]/(2u) gives the 
maximum of P in (|24|) with respect to r. 

Changing the variable u into t by u = cot t/v3, we find 
the entanglement E(cott/>/3) changes dynamically dur- 
ing the evolution as shown in Fig. || . It takes a value of 
at the initial average state |o) (t = n/6), attains its max- 
imum ~ 0.340 at the half-way state u = 1/3 (t = 7r/3), 
and finally goes back to at the target state \w) (t = 
7r/2). Note that Grover's algorithm in 2-qubits case uses 
entanglement at most ~ 0.340 though the available max- 
imal entanglement is n/2. This implies that, for the half- 
way state (t = 7r/3), there is a closer state on the quadric 
than either \a) or \w) whose distance from the half-way 
state is 2(7r/6) > 0.34. However, as seen in the following, 
Grover's algorithm comes to use the entanglement max- 
imally when the number of the qubit n becomes larger. 




t (rad) 

0.6 0.8 1 1.2 1.4 

FIG. 2. Entanglement E(cott/V3) in (M) during Grover's 
evolution for 2-qubits case (solid curve) is drawn, com- 
pared with the approximate estimate of the entanglement 
E 2 (cott/V3) in Q (dashed curve). Note that the approx- 
imate (dashed) curve agrees well with the exact (solid) one 
except near the half-way state (t = 7r/3). 



Now it is straightforward to calculate the entangle- 
ment for the general n-qubits case. According to ([20]), 
the Grover's state in the n-qubits case is given by 
(u, ■ ■ ■ ,u, 1), where u(:= cot(k+^)Q / N —I) ranges from 
to 1. The states \<f>), candidates closest to the state 
on the submanifold (C7 ,1 ) x ™, are "coherent" states (cf. 
(p3|)) parametrized as 




with v € C such that < |u| < 1. Likewise, by use of 
V = re %x with < r < 1, we have 



E n (u) = 2 arccos . /max P n {u, v) = 2 arccos 




2r cosx+ l) n + 2u(l - u) YT m =o (l) rm cos ™>x+0-- u) 2 
((iV-l)u 2 + l)(r 2 + l)™ ' 

(27) 



For a fixed r, the maximum in ( P7|) is attained at the 
phase x = (mod 27r) for the n-qubits case also. From 
■g-P n (u, r, X = 0) = 0, we have an extremum condition; 

u=- — ^ . (28) 

(1 + r)"- 1 (l - r) +r K ' 

It is hard to solve analytically the extremum condition 
( p8| ) for r so as to seek the solution that gives the maxi- 
mum of P n . However, we readily find, as seen in Fig. 
that ( pSj) increases monotonically with r for n < 6, on 
the other hand, it has a relative maximum and a relative 



minimum for n > 7, i.e. for almost all n. When u has 
one-to-one correspondence with r in Fig. ||, it soon be- 
comes the maximum condition of P n . In contrast, when 
u has one-to-three correspondence to r, the point, among 
the three, that is included in the solid parts of the curve 
in Fig. ^| indeed gives the maximum condition. Now it 
should be noted that because entanglement is symmetric 
for the half-way state as in the 2-qubits case, all we need 
to consider is one half of the whole dynamics, e.g. the sec- 
ond half here. By reparametrizing uasti = cot t/y/N—1, 
the second half, given by t G [(tt + 9)/4, 7r/2], corresponds 
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toue [cot[(7r + 0)/4]/\^V-T,O] 
the second half is almost u e [1/y/N 



Thus when N is large, 
111) so that 



1,0] by 

it can be treated as the realm of u < 1 and r <C 1. 
Taking first order of r in (P8|), we obtain an approximate 
maximum condition for r, u <C 1: u ~ + (n — l)r], or 



1 



(29) 



Though (|29j) becomes a better approximation for larger 
N, it seems to remain valid for small N because even in 
the 2-qubits (worst approximation) case , the deviation 
from the exact result is limited near the half-way state 
and is small (see Fig. |J). 



the entanglement of the n-qubits case: 



E n (u) ~ 2 arccos , 



[u(r M + 1)" + (1 - u)} 2 



[(N-l)u* + l]{r* M 



1)"' 



(30) 



drawn in Fig. || with u = cot t/ \/N — 1. We find in Fig. |] 
that E n {cott/y/N-l) almost converges to an enveloping 
triangle of E = -2\t - § | + § at n ~ 15. This sug- 
gests two points: first, entanglement is maximally used 
for large n. Second, the closest separable state during 
Grover's algorithm is either the initial average state \a) 
or the target state |u>), which implies the submanifold of 
completely separable states is sparse in the large n-qubits 
state space. 




Mhalf 



1 r 

FIG. 3. Schematic pictures of the extremum condition ( J2q ) 

for n < 6, n = nc := 4 + 2v2, and n > 7 are drawn. 
Mhalf := cot [(it + 9)/4]/y/N— 1 corresponds to the half-way 
state. When n > 7(> nc), the solid parts of the curve give 
the condition for the maximum of P n . 



IV. GEOMETRIC CONSTRUCTION OF 
OPTIMAL QUANTUM SEARCH 



In Sees. H and |l| we found that Grover's algorithm 
is a horizontal lift of a geodesic lying away from the sub- 
manifold of the separable states in CP, which can be 
interpreted as the geometric necessary condition for the 
optimal quantum search. In this section, let us consider, 
on the contrary, whether all the geodesies toward the tar- 
get state \w) become the optimal quantum search. That 
is to say, we discuss the geometric sufficient condition for 
the optimal quantum search, from which the bound of 
the computational time is also derived naturally. 



A. Geometric strategy by means of geodesies 



1 . 4 



1.2 




. 8 
. 6 
. 4 

t (rad) 

0.25 0.5 0.75 1 1.25 
FIG. 4. Each entanglement E n {cott/y/W^l) in @ for 
the n = 2, 3, • • •, 50-qubits case is drawn from the bottom to 
the top. Apparent singularities in the half-way states are just 
due to taking a mirror image of the approximate calculations 
( po| ) for the second half of the dynamics. While the true 
curves for small n should be smooth near the half-way states 
(cf. Fig. ^), an intrinsic singularity appears as the peak of an 
enveloping triangle E = — 2\t — ^| + ^ when n goes to infinity. 



Let us consider a set of all the geodesies through the 
target state \w) in CP^ -1 . As seen in Sec. |jj its hori- 
zontal geodesic in Ti. is just a real 2-dimensional rotation 
on the plane spanned by \w) and some arbitrary state \y) 
(denoted for brevity as the w-y plane). We can restrict 
\y) such that q (:— (w\y) € K) ranges from to 1, by 
choosing the preferable overall phase of \y) for each ray 
n(|y)). \w) and \y) are said to be "in phase" in terms 



of the Pancharatnan connection |14|. By a consequence 
of an elementary theorem of real Euclidean geometry, a 
2-dimensional rotation on the w-y plane is constructed 
by two successive reflections, 

Uy, := -Iyl w , (31) 

where I y i := 1 — 2|y')(y'|, I w := 1 — 2|w) (w| denote a 
reflection for the line orthogonal to \y'), \w) in the w-y 
plane, respectively. We take an overall —1 in (|l|) for 
convenience, which simply means that —I y i = I y > ± 
Though in general, \y') can be any state on the w-y 
plane, we put \y') = \y) without loss of generality be- 
cause \y) itself is any state. By using r\ £ [0,ir] s.t. 
(w\y), ( |3l| ) is represented by, 



Substituting tm of (E9h and \ = into (|27|), we have 



U V (V) 



-Iylu 



cos?y 
sinn 



- sm?7 
cos?7 



(32) 



G 



in the basis of \r')(:= (\y) — q\w))/^/l — q 2 ), orthogonal 
to |ui), and \w). Two remarks are in order: first, the 
angle of the rotation in (|3|), which corresponds to the 
speed of a single query, is determined just by 77, (or q). 
This means that the speed is faster for larger 77, (or q). 
Second, the direction of the rotation in (^ ) is determined 
by the order of I w and I y . Alternate applications of /,„ 
and I y cause successive rotations in the same direction, 
as can be seen in Fig. g. 

Thus the candidate for the algorithm that gives the 
optimal quantum search toward the target \w) is con- 
structed in terms of the geodesies as 



\^(k))=U Vk (r, k )---U y2 ( V2 )U yi ( m )\y } 



(33) 



such that |j/o), |yi), ■ • • , \yk) must lie on the same 2- 
dimcnsional plane including \w) with 770 < f]i < ■■• < ifk, 
where sin ^ := ( w \llj)- We find, however, only the case 
°f \Vo) = \yi) = ■■■ = \y k ) (i-c T70 = r?i = •••=%) is 
possible. 

The reason is the following. Suppose the algorithm 
begins from a fixed \yo), then = U yi (rji)\y ) is de- 

termined by selecting |yi) on the w-yo plane. However, 
because we never know (only the oracle knows) which 
is the target \w), the only state we are able to utilize 
on the w-yo plane is |yo)- So \yi) = \yo)- Then to get 
1^(2)) by the choice of £^(772), it might seem possible 
to apply |^(1)) as well as |yo) to \y 2 } on the w-yo plane. 
Yet, we must call another oracle as a subroutine to take 
advantage of |^(1)) in case no measurements are done 
during the computation. This case is just the situation 
where Zalka || showed the optimality of Grover's algo- 
rithm. We restrict our attention here to an algorithm 
including no subroutines which require another oracle, 
because, if needed, we can always embed no-subroutine 
algorithm into a certain larger algorithm as a subrou- 
tine |2E|j . Hence all we can do is \y 2 ) = \yo), again. In 
the same way, wc finally obtain \yo) = ■■■ = \y k ) (he. 
770 = ■■• = rjk) which can be denoted simply by \y) and 77, 
respectively. 




FIG. 5. Each horizontal geodesic toward the target \w) con- 
sists of 2-dimensional successive rotations in the plane includ- 
ing \w) . 



Accordingly, (|33|) turns out to be an extension of 
Grover's algorithm where the average state |o) is replaced 
with the arbitrary state \y) pt| . That is, our algorithm 
is written as 



\v) 



cos§|r') 
sin 



|*(fc)) := U y (r,) k \y) 



cos(k+^)r]\r'} 
sin(fc+|)r7|w) 



(34) 



The speed of the algorithm (Q) , considered as the trav- 
eling (Fubini-Study) distance of a single query along the 
geodesic, is given by 

As 

v (k) := -— = 2arccos|(*(fc+l)|*(fc))| = 2?y = 4arcsinq, 
Ak 

(35) 



where, in the third equality, we use (^(fc)^^)^^)) = 
cos 77 by ([52|) and (|34]). This corresponds to the Anandan- 
Aharonov relation |13[|: = 2AH/H, which means that 
the speed in CV is determined by the energy uncertainty 
AH. Note that V(k) in (35) is constant, independent of 



k, through the algorithm so that it depends only on 77 (or 
q). The total traveling distance is naturally thought to 
be the statistical distance [^7| between the initial state 
n(|y)) and the goal state n(|w)): 



7r — X] = tt — 2 arcsing. 



(36) 



Consequently, we obtain the time required to reach the 
target \w): 



w : ~ V 



7r — 77 7r — 2 arcsin q 



2'/ 



4 arcsin q 



(37) 



As seen in Fig. |], T w is shorter for larger q. We als o find 
that T w - 7r/(4g) for q < 1, while T w - ^2(1 - q)/ir for 
«~ 1. 



T 

1 11; 








'4-' q 



FIG. 6. The computational time T w in (|^). 
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B. Bound for the computational time 

We may ask where the bound of the computational 
time T w comes from. Remember that we want to extract 
\w) with probability 1 in an optimal computational time 
for the worst case evaluation. However, because we do 
not know a priori (only the oracle knows) which is the 
target \w), we have to select \y) independently of \w). 

When \y) is selected as \y) = J2x=o z x\ x ) (J2x=o \ z ^\ 2 = 
1) in the computational basis \x) (9 \w)), its smallest 
overlap q s ('-= \z s \) gives the computational time T w for 
the worst case. We find 

7r — 2arcsin<7, ir — 2arcsin(l/v7V) it j— 

— : — > wv^/ _ -VN, 38) 

4arcsmg s 4arcsin(l/V-/V) 4 



where because of q s < 1/vJV, the equality in (|38|) is 
attained for the \y) such that q s = 1/y/W, i.e., all 
Ix '■= \ z x\ = l/V^V- This implies that the "mixedness" of 
the searching state space (in part) bounds the efficiency 
of quantum search as Bose el al. 28 mentioned. Thus, 
regardless of which is the target, the optimal computa- 



tional time is js/N. This result of course coincides with 
Grover's result ||] , first proved optimal by Zalka || . 

It should be commented that there remains room for 
relative phases in \y). \y) = \a) in the original Grover's 
algorithm is only a choice. In general, any element of the 
Fourier basis, 

1 N ~ 1 

|P> : =^E e ^ P » (p = 0,l,...,iV-l), (39) 



can be taken as a \y) among the completely separable 
states. This implies that the quantum search takes ad- 
vantage of dual bases, \x) and \p), to run in the opti- 
mal computational time because its kernel takes the form 
U = -I P I X by d|). 



Grover's algorithm is exactly optimal regardless of n [EJ; 
(ii) the computational time is rather determined by the 
overlap q = (w\y) as seen in Sec. [fv|. It is significant that 
the algorithm consists of the shortest path by means of 
the geodesic; as a result it runs across entangled states 
away from the submanifold of the separable states. 

It is readily found that the multiple target case m is 
also characterized in completely the same manner. More- 
over, our geometric strategy would be useful to construct 
other efficient quantum algorithms, as some efficient clas- 
sical algorithms are widely known to be geodesies in their 
parameter spaces. Exploring the geometric viewpoint 
also seems appealing toward the realization of quantum 
computers. For instance, (i ) the holonomic approach 
to quantum computation p9], where loops by horizon- 
tal lifts of the path in CP construct the logic gates to 
compute quantum algorithms, is supposed to have built- 
in fault-tolerant features against local perturbations; and 
(ii) time optimal pulse sequences in NMR quantum com- 
puting pcfl, given by geodesies on certain coset spaces, 
would minimize the effect of relaxation and optimize the 
sensitivity of the experiments. 
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APPENDIX A: SEGRE EMBEDDING IN THE 
GENERAL CASE 



It is straightforward to extend the mapping / of the 
Segre embedding into the general case: 



V. CONCLUSIONS 

In this paper, we have shown two geometric char- 
acteristics of the quantum search: one is related to 
the geodesic in CP, the other related to entanglement. 
First, the geometric necessary and sufficient condition for 
the optimal quantum search is given by the horizontal 
geodesic joining the target \w) and a preferable selected 
initial state \y) such that it overlaps equally, up to rela- 
tive phases, with all the elements of the computational 
basis |x) (9 \w}). Second, Grover's quantum search uses 
entanglement for an arbitrary number of the qubits n, in 
particular almost maximally for large n. However, there 
seems to be no direct relationship between the amount 
of entanglement (how far the dynamics is away from the 
submanifold of separable states) and the optimal, i.e. 
shortest, computational time. This is because (i) the 
amount of entanglement is different for each n though 



/ : CP™ x CP" 1 ' — > CP( m+1 )( TO ' +1 ) -1 
((a , dm), (bo, fe m ')) ^ (a>obo,---,a b m >,aib , 

(Al) 

We find, in the same way as in the text, that the algebraic 
submanifold given by the image of / is the zero locus of 
all the homogeneous polynomials of degree 2: 

Q(i,j),(k,l) := z (m' + l)i+k z (m' + l)j+l ~ z {m' + l)i+l z (m' + l) j + k , 

(A2) 

where < i < j < m, < k < I < m! . Hence we have 

f(CV m x CP" 1 ' ) = {Q (iMk ,i) = 0}, (A3) 

where a set of quadratic constraints {Q(i,jY(fc,n = 0} con- 
sists of m(m + l)m'(m' + l)/4 simultaneous equations. 
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APPENDIX B: CALCULATION OF 
ENTANGLEMENT BY PARTIAL ENTROPY IN 
2-QUBITS CASE 



In this appendix, the entanglement of the 2-qubits case 
is calculated in terms of the partial entropy, so as to be 
compared with the results in Sec. III. For the Grover's 
state \ip) = (u|00)+u|01)+u|10} + |lT))/v / 3u 2 +l in @, 
we obtain a reduced density matrix by tracing out e.g. 
the second qubit; 

2u 2 
u(u 



p rcd := tr 2n d (IV'XV'I) 



1 



3u 2 + 1 



1 



u(u 

„,2 



(Bl) 

We calculate entanglement as the partial entropy of the 
reduced density matrix (Bl): 



e(u) := -tr(p rcd logp rcd ) = -A+ logA + - A_ logA_, 

(B2) 

where A+, the eigenvalues of the reduced density matrix 
p r ed in (pil), are given by 

i± vi-c(u) 2 



A+ : = 



C{u) := 2 



det- 



V3u r 



( u u\ 


2u(l - u) 


[u 1) 


3u 2 + 1 



(B3) 
(B4) 



As C(u) in (B4), called "concurrence" [|23j, goes from 
to 1, the entanglement e(u) in (B2) increases monotoni- 
cally from to 1. So C(u) as well as e(u) can be regarded 
as a measure of entanglement. 

The partial entropy e(cotf/-\/3) as well as the con- 
currence C(cott/v / 3), compared with our entanglement 
E(cott/y/3) in ( ^5|) , arc drawn in Fig. 0. We note 
two points in Fig. |7j: First, all the three entanglement 
measures are convex upward and are maximized at the 
half-way state t = tt/3 (u = 1/3). Second, surpris- 
ingly enough, our entanglement measure E almost co- 
incides with the concurrence C, characterized also as 
C = 2\z z 3 — ziz 2 \ — 2\Q\ (in (|l4|)) with normalized ho- 
mogeneous coordinates Zj. It should be noted, however, 
that their normalizations are different, i.e. E is normal- 
ized to 7r/2, while C is normalized to 1. In summary, our 
calculation of entanglement by (|2^) is consistent with 
that by the partial entropy. 



^/coti\ /~v cot £ \ tti/ cot t 



V3 



V3 



0.4r 
0.35.; 

0.3- 
0.25- 

0.2- 
. 15 

. 1 
0.05- 



t (rad) 



FIG. 7. Three entanglement measures: (i) the partial en- 
tropy e(cot t/y/S) in (B2) (the lower green solid curve); (ii) the 
concurrence C(cot£/v3) in (B4) (the upper red solid curve); 
and (iii) the minimum Fubini-Study distance to separable 
states _E(coti/\/3) in (G3) (the upper blue dashed curve). 
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